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Abstract 

I summarize the status of neutrino oscillations that follow from current data, including the 
status of the small parameters a = Arn^Q^/ AmJ^^^^ and sin'^ ^13 characterizing the strength of CP 
violation in neutrino oscillations. I briefly discuss the impact of oscillation data on the prospects 
for probing the absolute scale of neutrino mass in neutrinoless double beta decay. I also comment 
on the theoretical origin of neutrino mass, mentioning recent attemps to explain current data from 
first principles. 
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I. INTRODUCTION 



Nowadays neutrino physics lies at the center of attention of the particle, nuclear and 
astrophysics communities. Here I summarize the determination of neutrino mass and mixing 
parameters from current neutrino experiments, given in detail in Ref. [1]. (Waiting for the 
veredictum of MiniBoone, we neglect the LSND data). The key concept in these neutrino 
oscillation studies is that of neutrino mixing, a characteristic feature of gauge theories of 
massive neutrinos. The most complete study of the structure of the neutrino mixing matrix 
was given in [2]. For the analysis of current neutrino oscillation data one assumes its simplest 
unitary, 3-by-3, and CP conserving form, as there is currently no sensitivity to CP violation. 
The interpretation of the data requires good calculations of solar and atmospheric neutrino 
fluxes [3, 4], neutrino cross sections and response functions, as well the inclusion of matter 
effects [5, 6] in the Sun and the Earth. The reader is referred to Ref. [1] for technical details 
and further references. For the discovery prospects of future neutrino oscillation studies 
see the Neutrino Oscillation Industry Web- Page [7]. Testing for the effect of leptonic CP 
phases and the absolute scale of neutrino mass constitutes the main upcoming challenge. 
Dirac CP phases will be probed in future oscillations studies, while Majorana phases will be 
tested in future searches for PPqi, (neutrinoless double beta decay). I also briefly discuss the 
robustness of the oscillation hypothesis itself vis a vis the presence of non-standard physics. 

Last, but no least, where does the neutrino mass come from? I briefly comment on 
two alternative views on the theoretical origin of neutrino mass, mentioning some of their 
possible signatures. 

II. SOLAR AND REACTOR DATA 

The solar neutrino data includes the rates of the chlorine experiment (2.56 ± 0.16 ± 
0.16 SNU), the gallium results of SAGE (66.9 SNU) and GALLEX/GNO (69.3 ± 

4.1 ± 3.6 SNU), as well as the 1496-day Super-K data (44 bins: 8 energy bins, 6 of which 
are further divided into 7 zenith angle bins). The SNO sample includes the data from 
the salt phase in the form of the neutral current (NC), charged current (CC) and elastic 
scattering (ES) fluxes, the 2002 spectral day/night data (17 energy bins for each day and 
night period) and the 391-day data. The analysis includes both statistical errors, and 
systematic uncertainties such as those of the eight solar neutrino fluxes. 
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Reactor anti-neutrinos from KamLAND are detected at the Kamiokande site by the 
process Pg + P e"*" + n, where the delayed coincidence of the prompt energy from the 
positron and a characteristic gamma from the neutron capture allows an efficient reduction 
of backgrounds. Most of the incident PgS come from nuclear plants at distances of 80 — 350 
km from the detector, far enough to probe large mixing angle (LMA) oscillations. To avoid 
large uncertainties associated with geo-neutrinos an energy cut at 2.6 MeV prompt energy 
is applied for the oscillation analysis. 

The first KamLAND data corresponding to a 162 ton-year exposure gave 54 anti-neutrino 
events in the final sample, after cuts, whereas 86.8 ± 5.6 events are predicted for no oscil- 
lations with 0.95 ± 0.99 background events. This is consistent with the no-disappearance 
hypothesis at less than 0.05% probabihty, giving the first evidence for the disappearance of 
reactor neutrinos before reaching the detector, and thus the first terrestrial confirmation of 
oscillations with Amg^^^. 

Additional KamLAND data with a larger fiducial volume of the detector were presented 
at Neutrino 2004, corresponding to an 766.3 ton- year exposure. In total 258 events have 
been observed, versus 356.2 ± 23.7 reactor neutrino events expected in the case of no dis- 
appearance and 7.5 ±1.3 background events. This leads to a confidence level of 99.995% 
for Pg disappearance. Moreover they obtain evidence for spectral distortion consistent with 
oscillations. 

A very convenient way to bin the latest KamLAND data is in terms of rather 
than Ep,-. Various systematic errors associated to the neutrino fluxes, backgrounds, reactor 
fuel composition and individual reactor powers, small matter effects, and improved z/g flux 
parameterization are included [1]. Assuming CPT invariance one can directly compare the 
information obtained from solar neutrino experiments with the KamLAND reactor results. 

The KamLAND data single out the LMA solution from the previous "zoo" of alterna- 
tives [8] . The stronger evidence for spectral distortion in these data also leads to improved 
AmgQL determination, substantially reducing the allowed region of oscillation parameters. 
More than just cornering the oscillation parameters, however, KamLAND has eliminated all 
previously viable non-oscillation solutions. From this point of view KamLAND has played 
a key role in the resolution of the solar neutrino problem. 
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III. ATMOSPHERIC AND ACCELERATOR DATA 



The first evidence for neutrino conversions was the zenith angle dependence of the //-hke 
atmospheric neutrino data from the Super-K experiment in 1998, an effect also seen in other 
atmospheric neutrino experiments. At the time there were equally good non-oscillation 
solutions, involving non-standard interactions [9]. Thanks to the accumulation of upgoing 
muon data, and the observation of the dip in the L/E distribution of the atmospheric 
survival probability, the signature for atmospheric neutrino oscillations has now become 
convincing. The data include Supcr-K charged-current atmospheric neutrino events, with 
the e-like and //-like data samples of sub- and multi-GeV contained events grouped into 10 
zenith- angle bins, with 5 angular bins of stopping muons and 10 through-going bins of up- 
going muons. We do not use Ur appearance, multi-ring // and neutral-current events, since an 
efficient Monte-Carlo simulation of these data would require further details of the Super-K 
experiment, in particular of the way the neutral-current signal is extracted from the data. 
As far as atmospheric neutrino fiuxes are concerned, we employ the latest three-dimensional 
calculations given in [4]. 

The disappearance of i/^'s over a long-baseline probing the same Am^ region relevant 
for atmospheric neutrinos is now available from the KEK to Kamioka (K2K) neutrino os- 
cillation experiment. Neutrinos produced by a 12 GeV proton beam from the KEK proton 
synchrotron consist of 98% muon neutrinos with a mean energy of 1.3 GeV. The beam is 
controlled by a near detector 300 m away from the proton target. Comparing these near de- 
tector data with the content of the beam observed by the Super-K detector at a distance 
of 250 km gives information on neutrino oscillations. 

The first phase (K2K-1 data sample, corresponding to 4.8 x 10^^ protons on target) gave 
56 events in Super-K, whereas 80.1^54 were expected for no oscillations. The second phase 
(K2K-II data, corresponding to 4.1 x 10^^ protons on target) gave 108 events in Super-K, 
to be compared with 150.9lio!o expected for no oscillations. Out of the 108 events 56 are 
so-called single-ring muon events. This data sample contains mainly muon events from the 
quasi-elastic scattering -|- p — > /i + n, and the reconstructed energy is closely related to 
the true neutrino energy. The K2K collaboration also finds that the observed spectrum is 
consistent with the one expected for no oscillation only at a probabifity of 0.11%, whereas 
the best fit oscillation hypothesis spectrum has a probability of 52%. 

One finds that the neutrino mass-squared difference inferred from the i/^ disappearance in 
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K2K agrees with atmospheric neutrino results, providing the first terrestrial confirmation of 
oscillations with Am^TM with accelerator neutrinos. Due to low statistics the current K2K 
data sample gives a rather weak constraint on the mixing angle. However, although the 
determination of sin^ ^atm is completely dominated by atmospheric data, K2K data already 

start constraining the allowed Am^^^^, values [1]. For example, there is a constraint on Am^^^, 
from below, which is important for future long-baseline experiments, since their sensitivities 
are drastically affected if Am^TM lies in the lower part of the 3a range indicated by current 
atmospheric data. 



IV. THREE-NEUTRINO OSCILLATION PARAMETERS 

Lepton mixing is a characteristic feature of gauge theories of massive neutrinos. The first 
systematic study of the effective form of the lepton mixing matrix was given in [2]. For 
3-neutrinos the simplest form of this matrix can be taken as 

K = a;23<^i3<^i2 (1) 

where each factor in the product of the a;'s is effectively 2x2, characterized by an angle 
and a CP phase. Two of the three angles are involved in solar and atmospheric oscillations, 
so we set 9i2 = 6'sol and 6*23 = 6'atm- The last angle in the three-neutrino leptonic mixing 
matrix is ^13, 

/ ci3 e^-^^sia \ 
u;i3= 10 

V -e-'^''si3 ci3 / 
for which only an upper bound currently exists. All three phases are physical [10]. Two of 
the phases are fundamental, and arise at the two-generation level, being associated to the 
Majorana nature of neutrinos. They show up only in lepton-number violating processes, like 
neutrinoless double beta decay, not in conventional neutrino oscillations [10, 11]. The other 
phase corresponds to the phase present in the quark sector (Dirac-phase) and exists only 
with three (or more) neutrinos. This phase affects neutrino oscillations. 

Such unitary form for the lepton mixing matrix holds in models where neutrino masses 
arise in the absence of right-handed neutrinos. To a good approximation, it also holds if 
neutrino masses are induced by a high-energy-scale seesaw mechanism (see below). 

In our analysis we follow the simplest unitary form in Eq. 1. Since current neu- 
trino oscillation experiments are insensitive to CP violation, we also neglect all phases in 
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Figure 1: Three-neutrino regions allowed by the world's neutrino oscillation data at 90%, 95%, 
99%, and 3a C.L. for 2 d.o.f. In top panels is minimized wrt undisplayed parameters. 

the analysis. In this approximation oscillations depend on the three mixing parameters 
sin^ 9i2, sin^ 623, sin^ ^13 and on the two mass-squared splittings ArrigQ^ = ^"^21 = ^2 " ^1 
and Am^^j^j = Amg^^ = mg — ml characterizing solar and atmospheric neutrinos. The hier- 
archy ArrigQL ^ ^"^ATM implies that one can set Am^QL = 0, to a good approximation, in 
the analysis of atmospheric and K2K data. Similarly, one can set Am^^j^^ to infinity in the 
analysis of solar and KamLAND data. Apart from the data already mentioned, the analysis 
also includes the constraints from "negative" reactor experiments, CHOOZ and Palo Verde. 

The three-neutrino oscillation parameters that follow from the global oscillation analysis 
in Ref. [1] are summarized in Fig. 1 and in Table I. In the upper panels of the figure the Ax^ is 
shown as a function of the parameters sin^ 612, sin^ 623, sin^ ^13, Am2i, Amg^, minimized with 
respect to the undisplayed parameters. The lower panels show two-dimensional projections 
of the allowed regions in the five-dimensional parameter space. The best fit values and the 
allowed 3cr ranges of the oscillation parameters from the global data are summarized in 
Table I. This table gives the current status of neutrino oscillation parameters. 
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parameter 


best fit 


3a range 




e V J 


7.9 


7.1-8.9 




eV2] 


2.2 


1.4-3.3 


sin^ 6*12 




0.30 


0.24-0.40 


sin^ 6*23 




0.50 


0.34-0.68 


sin^ 6*13 




0.000 


< 0.043 



Table I: Neutrino oscillation parameters [1]. 




Figure 2: Determination of a = AnigQ^/ Am'^^j^ and bound on sin^ ^13 from current data, from [1]. 

Note that in a three-neutrino scheme CP violation disappears when two neutrinos become 
degenerate [2] or when one angle vanihes [12]. As a result its effects involve both the small 
mass hierarchy parameter a = Am^Qj^/Am^^^,, as well as the small mixing angle ^13. The left 
panel in Fig. 2 gives the parameter a, namely the ratio of solar over atmospheric splittings, 
as determined from the global analysis. The right panel in Fig. 2 gives Ax^ as a function 
of sin^ ^13 for different data samples. One finds that the KamLAND-2004 data have a 
surprisingly strong impact on this bound. Before KamLAND-2004 the overall bound on 
sin^ 6*13 was dominated by the CHOOZ reactor experiment, together with the determination 
of Amg^ from atmospheric data. 

In Fig. 3 we show the upper bound on sin^ ^13 as a function of Am^^^j from CHOOZ 
data alone compared to the bound from an analysis including solar and reactor neutrino 
data. One sees that, although for larger Am^^^,, values the bound on sin^ 613 is dominated 
by CHOOZ, this bound deteriorates quickly as Am^^^^ decreases (see Fig. 3), so that for 
Am?„., < 2 X lO-^eV^ the solar and KamLAND data become relevant. 
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sin^e,3 



Figure 3: Upper bound on sin^ ^13 (1 d.o.f.) from solar and reactor data versus Am^^^. Dashed 
(solid) curves correspond to 90% {3a) C.L. bounds, thick curves include KamLAND-2004 data, 
thin ones do not. Light (dark) regions are excluded by CHOOZ at 90% {3a) C.L. Hatched regions 
are excluded by ATM + K2K at 3a, horizontal line corresponds to current Am^^j^ best fit value. 

In summary the overall improvement is especially important for lower Am^^^ values. The 
new overall global bound on sin^ 613 is 0.043 at 3cr for 1 d.o.f. Such an improved sin^ ^13 bound 
follows mainly from the strong spectral distortion found in the 2004 KamLAND sample. 

Future long baseline reactor and accelerator neutrino oscillation searches [13], as well as 
studies of the day/night effect in large water Cerenkov solar neutrino experiments such as 
UNO or Hyper-K [14] could bring more information on sin^ ^13 [15]. With neutrino physics 
entering the precision age it is necessary to scrutize also the validity of the unitary approxi- 
mation of the lepton mixing matrix in future experiments, given its theoretical fragility [2] . 
Indeed, any model where neutrino masses follow "a-la-seesaw" gives corrections to this ap- 
proximation, which may be sizeable in some cases. 

V. ABSOLUTE NEUTRINO MASS SCALE 

Neutrino oscillation data are insensitive to the absolute scale of neutrino masses and 
also to the fundamental issue of whether neutrinos are Dirac or Majorana particles [10, 11]. 
On general grounds neutrino masses are expected to be Majorana [2], a fact that may 
explain their relative smallness with respect to other fermion masses. The significance of 
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the neutrinoless double beta decay stems from the fact that, in a gauge theory, irrespective 
of the mechanism that induces (5(5^^ , it necessarily implies a Major ana neutrino mass [16], 
as illustrated in Fig. 4. Hence the importance of searching for neutrinoless double beta 
decay [17]. Quantitative implications of the "black-box" argument are model-dependent, 
but the theorem itself holds in any "natural" gauge theory. 

Now that oscillations are experimentally confirmed we know that l3j3Qu must be induced 
by the exchange of light Majorana neutrinos, the so-called "mass-mechanism". The corre- 
sponding amplitude is sensitive both to the absolute scale of neutrino mass as well as the 
two Majorana CP phases that characterize the minimal 3- neutrino mixing matrix [2], none 
of which can be probed in oscillations. 

Fig. 5 shows the estimated average mass parameter characterizing the neutrino exchange 
contribution to jSlSoy versus the lightest neutrino mass. The calculation takes into account 
the current neutrino oscillation parameters in [1] and the nuclear matrix elements of [18] 
and compares with experimental sensitivities. The upper (lower) panel corresponds to the 
cases of normal (inverted) neutrino mass spectra. In these plots the "diagonals" correspond 
to the case of quasi-degenerate neutrinos [19, 20], which give the largest /3/3oi/ amplitude. In 
the normal hierarchy case there is in general no lower bound on the (5(5qv rate since there 
can be a destructive interference amongst the neutrino amplitudes. In contrast, the inverted 
neutrino mass hierarchy implies a "lower" bound for the /3/3oz/ amplitude. A normal hierarchy 
model with no lower bound on PPq^ is given in Ref. [21]. 

Future experiments like GERDA, SuperNEMO, CUORE, COBRA and others will extend 
the sensitivity and provide an independent check of the Heidelberg-Moscow claim [22, 23]. 
More information on the absolute scale of neutrino mass will also come from future beta 
decays searches (KATRIN) [24] as well as cosmology [25]. 




V V 



Figure 4: Neutrinoless double beta decay and Majorana mass are theoretically equivalent [16]. 
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Figure 5: Neutrinoless double beta decay amplitude versus current oscillation data. 

VI. THE ORIGIN OF NEUTRINO MASS 

It is well-known that the effective dimension- five operator ^£00 where the SU{2) f/(l) 
Higgs doublet and £ is a lepton doublet induces neutrino masses once the electroweak sym- 
metry breaks down through a nonzero vacuum expectation value (0) [26]. Nothing is known 
from first principles about the mechanism giving rise to this operator, its associated mass 
scale or flavour structure. The resulting Majorana neutrino masses are therefore unpredicted 
in general. However the very fact that Majorana neutrino masses violate lepton number may 
explain, irrespective of the underlying physics, why neutrinos are much lighter than the other 
fermions. 

One possibility is that the dimension-five operator is suppressed by a large scale in the 
denominator (top-bottom scenario). Alternatively, it may be suppressed by a small scale in 
the numerator (bottom-up scenario). Both scenarios are viable and can be made natural, 
the first being closer to the idea of unification. 
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Figure 6: Two types of seesaw mechanism. 

The most well-studied realization of the top-bottom scenario is the seesaw mechanism, 
which induces small neutrino masses from the exchange of heavy states that might come 
from unification. Small neutrino masses are induced either by heavy SU{2) ® f/(l) singlet 
"right-handed" neutrino exchange (type I) or heavy scalar bosons exchange (type II), in a 
nomenclature opposite from the one given in [2]. The effective triplet seesaw term has a 
flavor structure different from the type-I term, contributing to the lack of predictivity of 
general seesaw schemes, where both co-exist. Predicitvity within the seesaw approach may 
be obtained by appealling to extra symmetries, as given, for example, in [20]. The model 
predicts maximal 623 = vr/4, 613 = 0, and 6^12 = 0{1), though unpredicted. Moreover, if CP 
is violated ^13 becomes arbitrary but the Dirac CP violation phase is maximal [27]. The 
model leads to a variety of phenomenological implications. For example, it gives a lower 
bound on the absolute neutrino mass rrij, > 0.3 eV. It also requires a light slepton below 
200 GeV, and gives large rates for flavour violating processes. A survey of related models is 
given in Ref. [28]. 

Amongst "bottom-up" models we mention those where neutrino masses are given as ra- 
diative corrections [29] and models where low energy supersymmetry is the origin of neutrino 
mass [30]. The latter are based on the idea that R parity spontaneously breaks [31], leading 
to a very simple effective bilinear R parity violation model [32]. In this case the neutrino 
mass spectrum typically follows a normal hierarchy, with the atmospheric scale generated 
at the tree level and the solar scale radiatively "calculable" [33]. In order to reproduce the 
masses indicated by current data, typically the lightest supersymmetric particle decays in- 
side the detector. More strikingly, its decay properties correlate with neutrino mixing angles. 
For example, if the LSP is the lightest neutralino, it is expexted to decay 50/50 to muons 
and taus, since the observed atmospheric angle is close to 7r/4 [30, 33]. This constitutes a 
characteristic feature of the proposal that supersymmetry is the origin of neutrino mass and 
opens the tantalizing possibility of testing neutrino mixing at high energy accelerators, like 
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Figure 7: Flavour-changing effective operator for non-standard neutrino interaction, 
the "Large Hadron Collider" (LHC) and the "International Linear Collider" (ILC). 

VII. ROBUSTNESS OF THE OSCILLATION INTERPRETATION 

The general effective model-independent description of the seesaw at low-energies is char- 
acterized by (n, m), n being the number of SU{2) ® U{1) isodoublet and m the number of 
SU{2) ® U{1) isosinglet leptons [2]. This leads to a very rich and complex structure of the 
charged current lepton mixing matrix (non- unitary) and non-diagonal neutral current [2]. 
For example, the (3, 3) seesaw model has 12 mixing angles and 12 CP phases (both Dirac 
and Majorana-type) characterizing its full 3x6 lepton mixing matrix [2]. 

The nontrivial structure of charged and neutral current weak interactions is a general 
feature of seesaw models [2] and leads to dimension-6 terms non-standard neutrino interac- 
tions (NSI), as illustrated in Fig. 7. Such sub-weak strength sGp operators can be of two 
types: flavour- changing (FC) and non-universal (NU). In inverse seesaw-type models [34] the 
non-unitary piece of the lepton mixing matrix can be sizeable and hence the induced non- 
standard interactions may be phenomenologically important [35]. "Large" NSI strengths 
may also be induced by scalar boson exchanges in models with radiatively induced neutrino 
masses [29], and in supersymmetric unified models [36]. 

Non-standard physics may in principle affect neutrino propagation properties and detec- 
tion cross sections [37]. In their presence, the Hamiltonian describing atmospheric neutrino 
propagation has, in addition to the standard oscillation part, another term -ffNsi 



is the forward scattering amplitude for the FC process y^, + f i^t + f and y/^GpNfe' 





Here +(— ) holds for neutrinos (anti-neutrinos) and e and e' parameterize the NSI: \/2GFNfe 
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represents the difference between i^n + f and Ur + f elastic forward scattering. Here Nf is the 
number density of the fermion / along the neutrino path. In the 2-neutrino approximation, 
the determination of atmospheric neutrino parameters Am^^^ and sin^ ^atm was shown to 
be practically unaffected by the presence of NSI on down-type quarks (/ = d) [38]. Future 
neutrino factories will substantially improve this bound [39]. 

In contrast, the oscillation interpretation of solar neutrino data is more "fragile" with 
respect to the presence of non-standard interactions in the e — r sector [40]. On the other 
hand, it has been shown [41] that, even a small residual non-standard interaction of neutrinos 
in the e — r channel leads to a drastic loss in sensitivity in the ^13 determination at a neutrino 
factory. It is therefore important to improve the sensitivities on NSI, another window of 
oportunity for neutrino physics in the precision age. 

Work supported by Spanish grants FPA2005-01269 and BFM2002-00345 and by the EC 
RTN network MRTN-CT-2004-503369. 
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